Phase Diagram of the Extended Hubbard Model with Pair Hopping Interaction 
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A one-dimensional model of interacting electrons with on-site U, nearest-neighbor V , and pair- 
hopping interaction W is studied at half-filling using the continuum limit held theory approach. 
The ground state phase diagram is obtained for a wide range of coupling constants. In addition 
to the insulating spin- and charge-density wave phases for large U and V, respectively, we identify 
bond-located ordered phases corresponding to an enhanced Peierls instability in the system for 
W < 0, \U — 2V\ < \2W\, and to a staggered magnetization located on bonds between sites for 
W > 0, \U — 2V| < W . The general ground state phase diagram including insulating, metallic, and 
superconducting phases is discussed. A transition to the ^-superconducting phase at \U — 2V\ <C 
2t < W is briefly discussed. 

PACS numbers: 71.27.+a- Strongly correlated electron systems; heavy fermions- 71.10.Hf- Non- 
Fermi-liquid ground states, electron phase diagram and phase transitions in model systems- 71.10.Fd 
Lattice fermion models 



I. INTRODUCTION 

The one dimensional (ID) extended Hubbard model 
with nearest-neighbor repulsion V, in addition to the 
on-site repulsion U (hereafter U — V model) has been 
extensively studied during the last two decades as an im- 
portant theoretical test-bed for studying low-dimensional 
strongly correlated electron systems with rich phase 
structures. Considerable attention has been focused on 
stuying of the ground state (GS) phase diagram of the 
U — V model at half-filling, using analytical studies and 
numerical simulations jl]-[l5| The sketch of the phase di- 
agram consists of a Mott insulating phase (U > 2\V\) 
with dominating spin-density wave correlations, an in- 
sulating long-range-ordered (LRO) charge-density- wave 
(CDW) phase (2V > U > 0), and metallic phases with 
dominating singlet (SS) and triplet (TS) superconducting 
correlations. In the physically most interesting region of 
repulsive interactions (U, V > 0), the weak-coupling per- 
turbative renormalization group studies show that 
there is a continuous phase transition between SDW and 
CDW along the line U = 2V. In the strong coupling 
limit (U,V » 1) the SDW-CDW transition is discon- 
tinuous (first order) and the phase boundary is slightly 
shifted away from the line U = 2V ||,|]|§. Estimates 
for the location of the tricritical point, where the nature 
of the transition changes, have ranged from U c ~ 1.5 to 
U c ~ 5 (and V c ~ U c /2). §^§Jl| Recently increased 
interest towards the U — V Hubbard model was triggered 
by Nakamura, |ic|-[r^| who found numerically that for 
small to intermediate values of U and V, the SDW and 
CDW phases are mediated by the bond-ordered charge- 
density-wave (BO-CDW) phase. The SDW-CDW transi- 
tion splits into two separate transitions: (i) a Kostcrlitz- 
Thouless spin gap transition from SDW to BO-CDW and 
(ii) a continuous transition from BO-CDW to CDW. 



An analogous sequence of phase transitions in the 
vicinity of the U = 2V line is the intrinsic feature 
of extended U—V Hubbard models with bond-charge 
coupling ||T^ , [r7f . The bond located ordering in these 
models is directly connected with the site- off- diagonal 
nature of the bond-charge coupling. Models of corre- 
lated electrons with bond-charge coupling have currently 
attracted a great interest as models showing uncon- 
ventional, "kinematical" mechanisms of superconducting 
correlations |l8| [59|. Among the models with correlated 
"kinematics" models with pair-hopping interaction are 
the subject of current studies p7| -|39||. In this paper 
we consider the ground state phase diagram of extended 
U — V Hubbard model supplemented with the pair hop- 
ping term. The Hamiltonian of the model is given by 

T~i = ] ( c n,tr c n+l,<j + C n+l,(T C n,^) ~ A* ^ ] c n,cr c n,a 

n,(T n,a 

+ 7jU ^2 Pn,aPn,-<7 + V ^ PnPn+1 



(1) 



where p n>a = c\ a c n ^, p n = J2a P™,ai and (c„, ff ) 
denotes the creation (annihilation) operator for an elec- 
tron with spin a at site n. In Eq. (lj), t and fi denote 
the hopping integral and the chemical potential respec- 
tively, with U being the on-site Coulomb-Hubbard re- 
pulsion and V the intersite interaction. W is the pair 
hopping interaction. 

It is notable that the U, V, and W terms could be ob- 
tained from the same general tight-binding Hamiltonian 
p0| by focusing on a selected term of the two-particle 
interaction. The sign of the Coulomb-driven coupling 
constants is typically repulsive U,V,W > 0, and usually 
W <C U,V. However, below we will treat these param- 
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eters as the effective (phenomenological) ones, assuming 
that they include all the possible renormalizations, and 
their values and signs could be arbitrary. 

Interest in models with pair-hopping coupling comes 
from the unusual mechanisms of Cooper pairing provided 
by this interaction. In the absence of the on-site and 
nearest- neighbor couplings (U = V = 0), the model Eq. 
(0) reduces to the Penson-Kolb (PK) model @. The 
PK model is possible the simplest model which captures 
the essential physics of an electron system showing the 
ry-superconductivity in the ground state. In the 77-paired 
state, the eigenstates of the correlated electrons are con- 
structed exclusively in terms of doublon (on-site singlet 
pair) creation operators 0|. We consider two different 
realizations of the 77-paired state, constructed in terms 
of zero size Cooper pairs with center- of -mass momentum 
equal to zero (770-pairing) and 7r (^-pairing) , respectively. 

In the case of an "attractive" (W < 0) pair-hopping 
interaction the PK model describes a continuous evo- 
lution of the usual BCS type superconducting state at 
|J7| , \ W\ <C t into a local pair ^-superconducting state 
at \W\/t ->■ 00 @. In the case of repulsive (W > 0) 
pair-hopping interaction, in contrast, the transition into 
the ^-paired state takes place at finite W c and is of first 
order (level-crossing type) |3^j3^,^] . 

In this paper we address the question, whether the 
pair-hopping coupling could lead to the superconduct- 
ing ordering in the physically most relevant region of 
parameters U, V ^> W > 0. In this communication we 
present the weak-coupling ground state phase diagram of 
the model Eq. ([[]). As we show in this paper, the "attrac- 
tive" (W > 0) pair-hopping coupling enlarges the region 
of coupling constant corresponding to the metallic phase 
with dominating SS and TS instabilities. However, in 
the repulsive sector of the phase diagram, along the line 
U = 2V > 0, only the insulating LRO BO-CDW phase, 
is realized at 1 1/ — 2V^| < \W\. In the case of a "repulsive" 
(W > 0) pair-hopping coupling, the BO-SDW phase cor- 
responding to a bond located staggered magnetization is 
together with the CDW the most divergent instability 
in the system. We also present quantitative arguments 
in favour of an additional phase transition at W ~ 4£ 
from the insulating BO-SDW to the ^-superconducting 
phase. 

The outline of the paper is as follows: In section II 
we present the continuum limit bosonized version of the 
model. In section III we discuss the weak coupling phase 
diagram. Section IV is devoted to a discussion of the 
ground state phase diagram and a summary. 

II. CONTINUUM LIMIT THEORY AND 
BOSONIZATION. 

In this section we derive the low-energy effective field 
theory of the lattice model Eq. (1) at half-filling. Con- 



sidering the weak-coupling case |Z7|, |V|,|W| <C t , we 
linearize the spectrum and pass to the continuum limit 
by use of the mapping 

% 1/2 c„, CT - i n R a (x) + {-i) n L a {x). (2) 

Here x = nag, ao is the lattice spacing, and R a (x) and 
Ltj(x) describe right-moving and left-moving particles, 
respectively. These fields can be bosonized in a standard 
waygi: ' 

R a (x) = {2na r 1/2 e l ^^^ x \ (3) 
L a (x) = (27ra r 1/2 e~ lvq "*^ (a;) , (4) 

where &r(l),<t(x) ar e the right (left) moving Bose fields. 
We define $ CT = &R, a + $l,<t and introduce linear com- 
binations, tfc = ($ T + $;)/V2 and tp s = ($ T - $|)/\/2, 
to describe the charge and spin degrees of freedom, re- 
spectively. Then, after a rescaling of fields and lengths, 
we rewrite the bosonized version of the Hamiltonian 
(|l|) in terms of two decoupled quantum SG theories, 
H = H C +H 8 , where 

H e(s) = fdx{^ [{d x ^ c{s) f + (d^ c , (s) ) 2 ] 

+ ^^cos(^8irK c{s) cp c (x))y (5) 

Here C ( S )( X ) are the dual counterparts of the fields 
</ ) c(s)( ;r ) : dx9 c (s) = II C ( S ) where n c ( s ) is the momentum 
conjugate to the field 4> c ( s ) ■ Here we have defined 

K c = (1 - g c )- 1 ' 2 ~ 1 + \g e , m e = -p-, (6) 
K s = (1 - g s )- 1 ' 2 si 1 + \g„, m s = ^, (7) 

v c (s) = v fK~^ are the velocities of the charge and spin 
excitations, Vf = 2tao(l — W/nt), and the small dimen- 
sionless coupling constants are given by 

g s = g x = (U-2V+2W)/nv F , (8) 
g c = -(U + 6V + 2W)/ttv f , (9) 
g u = (U - 2V - 2W)/ttv f . (10) 

The relation between K c (K s ), m c (rn s ), and g c (g s ), g u 
(gj_) is universal in the weak coupling limit. 

In obtaining (|^) the strongly irrelevant term ~ 
cos(v / 87r K c if c ) cos(\/8n K s ip s ) describing umklapp scat- 
tering processes with parallel spins was omitted. The 
mapping of the Hamiltonian (Q) onto the quantum the- 
ory of two independent charge and spin Bose fields, allows 
a study of the ground state phase diagram of the initial 
electron system using the far-infrared properties of the 
bosonic Hamiltonians (||) . Depending on the relation be- 
tween the bare coupling constants K and m the infrared 
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behavior of the quantum SG field exhibits two different 
regimes Q: 

For |m| < 2{K — 1) we are in the weak coupling regime; 
the effective mass M — > 0. The low energy (large dis- 
tance) behavior of the gapless charge (spin) excitations 
is described by a free scalar field. The corresponding- 
correlations show a power law decay 

( e W2*K * v {x) e -iV2wK*<e(x')} _ \ x _ x l\~K* ^ ^ 

^^•si^-Vh/x'JW) „ | x _ ^j" 1 /** , (12) 

and the only parameter controlling the infrared behav- 
ior in the gapless regime is the fixed-point value of the 
effective coupling constants K*( s y 

For \m\ > 2 (K — 1) the system scales to a strong cou- 
pling regime: Depending on the sign of the bare mass to, 
the effective mass M — > ±oo, which signals the crossover 
into a strong coupling regime and indicates the dynam- 
ical generation of a commensurability gap in the exci- 
tation spectrum. The field y c ( s ) gets ordered with the 
vacuum expectation values E3] 



, .{^m^ (to> ) 

\<P*'V ~\o (to < 0) 



(13) 



The ordering of these fields determines the symmetry 
properties of the possible ordered ground states of the 
fermionic system. 

Using Eqs. (|§|)-([l0|) and ([lj), one easily finds that there 
is a gap in the spin excitation spectrum {M s — > — oo) for 

U -IV + 2W <0. 

In this sector of coupling constants, the ip s field gets 
ordered with vacuum expectation value (tp s ) = 0. At 
U — 2V + 2W > the spin excitations are gapless and 
the low-energy properties = of the spin sector are de- 
scribed by the free Bose field system with the fixed-point 
value of the parameter K* = 1 . 
The charge sector is gapped for 



and for 



U > max{2V + 2W,-2\V\} 



U <2V + 2W but 2V + W > 0. 



with the fixed-point value of the parameter 
K* ~ 1 + y/2(U + 2 V) (W + 2 V) /irv F . 



(14) 



Especially important is the line U = 2V + 2W corre- 
sponding to the fixed-point line m c = 0, K c — 1 < 0. 
Here the infrared properties of the gapless charge sector 
are described by a free massless Bose field with the bare 
value of the Luttinger liquid parameter K c . 

To clarify the symmetry properties of the ground states 
of the system in different sectors we introduce the follow- 
ing set of order parameters describing the short wave- 
length fluctuations of the 

• s^e-located charge and spin density: 



Acw = (-l)"J>„, 



±SDW 



sin ( y/ 2nK c ip c ) cos ( \/2i^K a ip a ) 



~ cos(\/2iiK c (p c ) sin(^/2nKs^p s ) , 
• bond-located charge-density: 

&-BO-CDW = (-l)"^(4^C n+ i !ff + h.C.) 



(15) 



(16) 



~ cos(v/2tt K c (p c ) cos(y/2irK s (p s ) 
• The bond-located spin-density: 

Abo-sdw = (-1)" cr(cj lj0 .c n+ i !0 - + h.c.) 



(17) 



sm.(y/2n K c <Pc) sm{y/ 2n K s tp s ) . (18) 



In addition we consider two superconducting order pa- 
rameters corresponding to the 

• singlet and triplet superconductivity: 



Ass(x) = R\{x)L\(x) - r\(x)l\(x) 



exp(i\ —0 C ) cos(\/2tt K s ip s ), 
V A r 



A TS (a;) = R\(x)L\(x) + r\(x)l\{x) 



exp(iw— 6 C ) sm(y / 2wK s (p s 



(19) 



(20) 



In the former case M c — > — oo and the vacuum expecta- 
tion value of the charge field (ip a ) = 0, while in the latter 
case M c — » oo and (<p s ) = ^n/8K c . 

In the sectors of coupling constants corresponding to 
the gapless charge excitation spectrum the properties of 
the charge degrees of freedom are described by the free 
Bose field 



Vc 

2 



K, 



III. WEAK-COUPLING PHASE DIAGRAM 

With these results for the excitation spectrum and the 
behavior of the corresponding fields, Eqs. (|Tl|)-(|l3|), we 
now discuss the weak- coupling ground state phase dia- 
gram of the model (Q) . Below we will focus on the new 
phases appearing in the phase-diagram due to the effect 
of the pair-hopping coupling. The phase diagram consists 
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of 5 sectors (see Fig. 1 and Fig. 2). Sectors A, B, CI, C2 
are present in the phase diagram of the U — V Hubbard 
model @]. 
Sector A 

• U> max{2V + 2W,-2\V\}, 

corresponds to the ordinary Mott insulating phase: The 
charge excitation spectrum is gapped, the spin sector is 
gapless. The ordering of the field tp c with vacuum ex- 
pectation value (ip c ) = leads to a supression of the 
superconducting, CDW, and BO-SDW correlations. The 
SDW and Dimer correlations show a power-law decay at 
large distances 

(Asdw(x)Asdw(x')) ~ (Abo-cdw(x)Abo-cdw(x')} 
r^lx-x'f 1 . (21) 

Sector B 

• U < 2V - 2\W\ and 2V + W > 

corresponds to the long-range ordered CDW insulating 
phase. The charge and spin excitations are gapped. The 
fields <p c (s) S e t ordered with vacuum expectation values 

c , and 



and (tp c ) = 
(A C dw(x)A C dw(x')) 



constant. 



(22) 



Sector CI 



• U < min{2V - 2W; ~2V} and 2V + W < 

corresponds to the Singlet Superconducting (SS) phase. 
There a gap exists in the spin excitation spectrum and 
the spin field is ordered with (<^ s ) = 0. The charge ex- 
citation spectrum is gapless with the fixed point walue 
of the parameter K* > 1. The SDW, BO-SDW, and the 
TS instabilities are suppressed. The CDW, BO-SDW, 
and the SS instabilities show a power-law decay at large 
distancesr. However since K* > 1 the SS instability 



(A ss (x)A ss (x')) ~ \x-x'\ 

dominates in the ground state. 
Sector C2 



-l/K e 



(23) 



• -2\V\ < U < 2V - 2W and 2V + W < 

corresponds to the Luttinger liquid phase with dominat- 
ing superconducting instabilities. None of the conditions 
of charge and spin gap is satisfied here. In this sector, 
the system shows the properties of a Luttinger liquid with 
dominating superconducting instabilities TS and SS. The 
singlet superconducting and triplet superconducting cor- 
relations show the same power law decay at large dis- 
tances and the TS instability dominates because of the 
weak logarithmic corrections n. 

Finally we analyze the sectors describing the new 
phases. These new phases essentially appear along the 




FIG. 1. The ground state phase diagram of the ID 
U — V — W model for the case of a half-filled band and W < 0. 
Solid lines separate different phases: A. SDW + BOW - 
Mott insulating phase with an identical power-law decay of 
spin-density-wave and Peierls correlations. B. CDW (LRO)- 
long range ordered (LRO) charge density wave phase; CI. 
Singlet superconducting phase. C2. Metallic phase with dom- 
inating singlet and triplet superconducting correlations. D. 
LRO dimerized (Peierls) phase. 



SDW-CDW transition line U = 2V > of the U - V 
Hubbard model. In the weak-coupling limit, the transi- 
tion from the Mott insulating phase at U > 2V to the 
CDW insulator at U < 2V is mediated by the Luttinger 
liquid phase with gapless spin and charge excitations. 
AT U = 2V the Mott insulator charge gap closes and at 
U — 2V < the charge and the spin gap opens simulta- 
neously. In the very presence of the pair-hopping inter- 
action, the SDW-CDW transition splits into two transi- 
tions: Along the line U = 2V — 2W the spin gap opens, 
while at U = 2V + 2W the Mott insulator charge gap 
closes, and for U < 2V + 2W the CDW charge gap 
opens. In the case of an attractive pair- hopping inter- 
action W < (Fig. 1) the spin gap opens in the presence 
of a Mott insulator charge gap. Therefore, in sector D 



\U-2V\ < 2\W\ 



2V + W > 0, 



the charge and spin channels are gapped and both, charge 
and spin fields are ordered, (ip c ) = (ip a ) = 0. In this case 
the long-range ordered BOW phase 



(A 



BO-CDW 



(x)A 



BO-CDW 



(*')> 



constant 



(24) 



is realized in the ground state. 

In the case of a repulsive pair-hopping coupling W > 
(Fig. 2), the transition within the charge degrees of free- 
dom, takes place before the spin gap opens. Therefore, 
in sector Dl 
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FIG. 2. The ground state phase diagram of the ID 
U — V — W model for the case of a half-filled band and 
W > 0. The dot-dashed line marks a transition to the 
^-superconducting phase. 

• \U-2V\ < 2\W\ and U + 2V > 0, 

the generation of a gap in the charge excitation spectrum, 
accompanied by the ordering of the field ip c with vacuum 
expectation value (ip c ) = \J n/8K c , leads to a supression 
of the superconducting, SDW, and BO-CDW ordering. 
The CDW and BO-SDW correlations show a power-law 
decay at large distances 

(Acdw(x)Acdw(x')) ~ (Abo-sdw(x)Abosdw(x')) 

Hz-aT 1 . ( 25 ) 

Therefore, this sector of the phase diagram corresponds 
to the insulating phase with coexisting CDW and BO- 
SDW instabilities. 

Let us now discuss the ^-superconducting phase. In 
models with "kinematical" mechanisms of Cooper pair- 
ing, transition to an 77-paired phase is typically the finite- 
bandwidth phenomenon |l9| , p2| ^4j. In the case of pair- 
hopping interaction, the transition point is determined 
by the competition between the single-electron and dou- 
blon derealization energies. After the transition the con- 
tribution of the one-electron hopping term to the ground 
state energy almost vanishes and the ground state en- 
ergy is determined by the created strongly correlated 
two-particle ^-pair band |35|,|3{|. Simultanously, after 
the transition the spin gap opens in the system while the 
charge gap (at half-filling) closes (3lj . In the case of the 
PK model the transition point W C (U = V = 0) ~ 1.8* 
p5| , while in the case of the on-site Hubbard repul- 
sion, W C (V — 0) ~ 1.8t + all, where a is of the or- 
der of unity [[39). In both cases the insulating CDW 



+ (BO-SDW) phase is unstable toward transition to the 
^-superconducting state [ j36|]3^ , |39| ] . Due to the finite- 
bandwidth nature of the transition to a ^-paired state, it 
could not be consistently studied within the continuum- 
limit (infinite band) approach. Nevertheless, the exis- 
tence of a transition is clearly traced in the additive 
renormalization of the Fermi velocity (bandwidth) by the 
pair-hopping term vf = 2tao(l — W/nt). In the nar- 
row stripe along the frustration line \U — 2V\ <C W, 
the effects of the on-site and nearest-neighbor repulsion 
cancel each other. The dimensionless coupling constants 
controlling the spin degrees of freedom (||) are exactly 
the same as in the case of the PK model. Therefore 
we conclude that along the frustration line U = 2V 
an additional phase transition with increasing W from 
the BO-SDW to the ^-superconducting takes place with 
W c ~ W C (U = V = 0) ~ 2t. Numerical studies of this 
sector of the phase diagram are currently in progress and 
will be published elsewhere. 

IV. DISCUSSION AND SUMMARY 

To summarize, we have presented the weak-coupling 
ground state phase diagram for ID extended U —V Hub- 
bard with pair-hopping in the case of a half-filled band. 
We have shown that the model has a very rich phase di- 
agram which includes the singlet-superconducting phase, 
a metallic phase with dominating SS and TS instabil- 
ities and four different insulating phases corresponding 
to the Mott antiferromagnet, the CDW insulator, the 
bond-ordered CDW and the bond-ordered SDW phase. 
In addition, we argued for the existence of a phase transi- 
tion to the ^-superconducting phase within the narrow 
stripe at \U - 2V\ « 2t < W. 
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